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ABSTRACT 



Context. Circumstellar disks and envelopes of low-mass young stellar objects (YSOs) contain significant amounts of ice. Such icy 
material will evolve to become volatile components of planetary systems, such as comets in our solar system. 

Aims. To investigate the composition and evolution of circumstellar ice around low-mass YSOs, we observed ice absorption bands in 
the near infrared (NIR) towards eight YSOs ranging from class to class II, among which seven are associated with edge-on disks. 
(— I ' Methods. We performed slit-less spectroscopic observations using the grism mode of the Infrared Camera (IRC) on board AKARI, 

^ ^, which enables us to obtain full NIR spectra from 2.5 /jm to 5 yum, including the CO2 band and the blue wing of the H2O band, 

I . which are inaccessible from the ground. We developed procedures to carefully process the spectra of targets with nebulosity. The 

O ' spectra were fitted with polynomial baselines to derive the absorption spectra. The molecular absorption bands were then fitted with 

the laboratory database of ice absorption bands, considering the instrumental line profile and the spectral resolution of the grism 
C/5 , dispersion element. 

, Results. Towards the class O-I sources (LI527, IRC-LI041-2, and IRAS04302), absorption bands of HjO, CO2, CO, and XCN are 

clearly detected. Column density ratios of CO2 ice and CO ice relative to H2O ice are 2 1 - 28 % and 1 3 - 46 %, respectively. If XCN is 
OCN", its column density is as high as 2 - 6 % relative to H2O ice. The HDO ice feature at 4. 1 fjm is tentatively detected towards the 
class O-I sources and HV Tau. Non-detections of the CH-stretching mode features around 3.5 pm provide upper limits to the CH3OH 
abundance of 26 % (L1527) and 42 % (IRAS04302) relative to H2O. We tentatively detect OCS ice absorption towards IRC-L1041-2. 
Towards class O-I sources, the detected features should mostly originate in the cold envelope, while CO gas and OCN" could originate 
in the region close to the protostar, where there are warm temperatures and UV radiation. We detect H2O ice band towards ASR4I 
and 2MASSJ 1628 137-243 139, which ai'e edge-on class II disks. We also detect H2O ice and CO2 ice towards HV Tau, HK Tau, and 
UY Aur, and tentatively detect CO gas features towards HK Tau and UY Aur. 
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1 . Introduction heavily extincted by dust, and thus are faint. Secondly, it is un- 

. ^ , clear whether the ice absorption bands, if detected, originate in 

^ . In molecular clouds, protostellar envelopes, and protoplanetary the disk or other foreground components (i.e. ambient clouds). 

^ disks, significant amounts of oxygen and carbon are in the form The source CRBR 2422.8-3423 was the first edge-on object to- 

_C3 ^ of molecules in ice m antles, such as H 2O, CO , CO2, and CH3OH wards which detailed ice o bserva tions were performed (Thi et al. 

(e.g. Whittet 119931 Murakawa et al. 120001 Gibb et al. I2004I) . |2002b . Pontoppidan et al. ( l2005l l concluded, via detailed analy- 

These interstellar ices are formed by the adsorption of gas-phase sis and modeling of the object, that only a limited amount (< 20 

molecules onto grain surfaces and/or the g rain-s urface reactions %) of detected CO ice may exist in the disk, while up to 50 % 

of the adsorbed species (e.g. Aikawa et al. l2005b . of water and CO2 may originate in the disk. They also found 

Ices in circumstellar disks of low-mass YSOs are of spe- that the 6.85 yum band, which is tentatively attributed to NH+, 

cial interest as they contribute to the raw material of planetary has a prominent red wing. Since this wing is reproduced by the 

formation. Observations of disk ices, however, are not straight- thermal processing in the laboratory, it indicates that along the 

forward for the following two reasons, and thus remain Umited line of sight towards CRBR 2422.8-3423 there is warm ice in the 

in number Firstly, ices exist mostly near the midplane at r > disk. 

several AU. Observations of disk ices are basically restricted to Honda et al. (120091 ) succeeded in detecting water ice in a 

edge-on disks, where the ice bands absorb stellar light, scattered disk around a Herbig Ae star, which is not edge-on. They ob- 

stellar light, and/or thermal emission of warm dust in the disk served scattered light from a disk around HD 142527 at multiple 

inner radius (Pontoppidan et al i2005i) . The targets are naturally wavelength bands around ~ 3 fim. This method can be a power- 
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ful tool for investigating the spatial distribution of ices in disks. 
However, it would cuiTently be difficult to detect ice features 
other than the 3 //m water band via scattered light. 

In the present paper, we report near infrared (NIR) (2.5 - 5 
fim) observations of edge-on disks using the grism mode of the 
IRC on board AKARI. We selected seven edge-on young stellar 
objects (YSOs); three are considered to be in a transient state 
of class O-I, and four are class II (see §2 for more detailed de- 
scriptions). Although the absorption by the protostellar envelope 
may overwhelm the absorption by the disk in class O-I YSOs, 
envelope material along our line of sight (perpendicular to the 
rotation axis and outflow) is likely to end up in the midplane 
region of the forming disk, and ices with a sublimation temper- 
ature of > 50 K (i.e. CO2 and H2O) would survive the accretion 
shock onto the disk (Visser et al. 12009} . 

Ice observations towards low-mass YSOs have been per- 
formed using Spitzer Space Telescope ( SST) and various ground- 
based telescopes such as the Very Large Telescope (VLT) and 
Subaru. Pontoppidan et al. (I2003I I observed the 4.67 fim CO ice 
band, which is observable from the ground, towards 39 YSOs. 
Using SST, Furlan et al. (2008) and Pontoppidan et al. ( l2008l l 
observed the strong CO2 ice band at 15.3 //m towards class I 
sources. Observations with SST also detected absorption bands 
of CH3OH, H COOH H2CO, and CH4 (Boogert et al. l2008l 
Zasowski et al. 120091 Oberg et al. '2008). Since ground-based 
observations are restricted to atmospheric windows, and since 
SST is restricted to the mid infrared (5 - 36 //m), NIR observa- 
tions via AKARI are important to comprehend the full spectra of 
ices in circumstellar material. AKARI has a high enough sensi- 
tivity for the spectroscopic observation of low-mass YSOs; the 
sensitivity of the AKARI grism mode at 3 /jm is 120 /iJy (1 cr, 
lOmin), which is almost comparable to the K-band sensitivity 
of the spectroscopic observation at 8-m ground-based telescopes 
(VLT and Subaru). While the L-band sensitivity on the ground 
is ~ 5 mag or 100 times worse than at K-band, the sensitivity 
of the AKARI grism mode does not vary significantly at 2.5 - 4 
jum, and only gradually worsens at longer wavelengths, by up to 
a factor of 3-4 at 5 fim. 

2. Observations 

Observations were performed using the IRC (Onaka et al 120071 
Ohyama et al. 12007) on board the AKARI satellite (Murakami 
et al 200711 from September 2006 to March 2007. We used the 
Astronomical Observational Template 04 (AOT04), designed 
for spectroscopy. The AOT04 replaces the imaging filters with 
transmission-type dispersers on the filter wheels to take near- 
and mid-infrared spectra. The field of view (FOV) consists of a 
r X r aperture (covered by 40 x 40 pixels in the imaging mode) 
(Figure [T^) and a 10' x 10' aperture. The targets are set at the 
center of the 1' x 1' aperture (without a slit), where only the NIR 
channel is available. The resolving power of the IRC grism mode 
is 0.00965 //m/pixel, and the FWHM of the point spread func- 
tion (PSF) is ~ 2.9 pixels in the NIR. The spectral resolution 
is thus ~ 100, which is sufficient to detect the ice absorption 
bands. In each observation, the following data are taken in a pe- 
riod of 600 s: two dark frames, one image frame with the N3 
filter (2.7 - 3.8;um), and eight slitless spectral images. 

Table [T] lists our targets, and Figure [T] (a) shows the N3 im- 
ages of our targets. The stellar light is dispersed from right (short 
wavelengths) to left (long wavelengths). The position angle of 
the FOV is fixed; the dispersion direction is almost aligned with 
the ecliptic latitude. The nebulosity of our targets extends almost 
perpendicularly to the dispersion direction. As an example, the 



spectral images of L1527 and IRAS04302 are shown in Figured] 
(b). 

Both IRAS04302 and L1527 (IRAS 04368H-2557) are edge- 
on YSOs in a transient state from class to class I (Andre et al. 
20001 . The YSO L1527 is of special interest from a chemical 
point of view; it is one of two YSOs where high abundances of 
carbon-chain species are observed. Carbon chains are considered 
to be formed by warm carbon-chain chemistry (WCCC): gas- 
phase reactions of sublimated CH4 (Sakai et al. 12008! Aikawa 
et al. 2008). The large organic species detected towards several 
low-mass YSOs (e.g. IRAS 16293) are not, however, detected to- 
wards L1527 (Sakai et al. 2008). Since the large organic species 
are formed from CH3OH, while carbon chains are formed from 
CH4, it is important to investigate ice bands, including the C-H 
bands, towards LI 527, to determine whether its ice composition 
significantly differs from those of other low-mass YSOs. 

The YSO IRC-L1041-2 was initially in our target fist of 
background stars (i.e. field stars behind molecular clouds) 
(Noble et al. in prep). The N3 image shows a faint object with 
nebulosity, which is similar to LI 527. The 557" observations also 
detected a YSO with nebulosity (T. Bourke private communica- 
tion). The presence of a dark dust lane indicates that the YSO 
has a high inclination, close to edge-on. 

The YSOs 2MASS J1628137 and ASR41 are edge-on class 
II objects. 2MASS J1628137 is accompanied by a tenuous 
(4 X 10""* Mq) envelope (Grosso et al. 2003). Hodapp et al. 
(120041) performed detailed modeling of the scattering and radia- 
tive transfer of ASR41 and concluded that the K-band image can 
be well reproduced by a protoplanetary disk embedded in a low- 
density cloud. 

HK Tau B and HV Tau C are edge-on class II objects, in 
the HK Tau binary and HV Tau triplet systems, respectively. We 
initially aimed to obtain spectra of these edge-on objects, but 
the spectra were found to be contaminated and dominated by the 
bright primaries, because the binary separations are small (see 
§4.1 and §5). The ice absorption features towards HK Tau B 
and HV Tau C might still be observable, although the derivation 
of ice column densities is not straightforward. For a comparison 
with disks that are not edge-on, we also observed UY Aur, which 
is a binary system with a separation of 0.9" accompanied by a 
circumbinary disk. The inclination of the disk is ~ 42° (Hioki et 
al. l2007l l. 

3. Data reduction and analysis 

3.1. Reduction of the spectral image 

The reduction pipeline provided by the AKARI project team, 
IRC.SPECRED (Ohyama et al l20U7]) . is unsuitable for our tar- 
gets, because it automatically defines the area suiTounding se- 
lected sources as "sky", an area that is often contaminated by 
nebulosity in our observations. 

To derive spectra from the data, we developed our own pro- 
cedure, referring to the reduction method by Sakon et al. ( 120071 1. 
Firstly, we made a cosmic-ray-free dark for each target by tak- 
ing the median of three dark frames; two dark frames were taken 
during the observation of each target, and one additional dark 
frame was taken from the preceding or subsequent observation 
in the AKARI archive. This procedure was necessary, because 
some frames are heavily damaged by cosmic rays. The dark was 
subtracted from the N3 image and spectral images. 

The eight spectral images were slightly shifted in relation 
to each other by the attitude drift of the satellite. The task 
IRC-SPECRED measures the shift of each image relative to the 
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Name 


coordinates (12000)" 


Class 


note 


L1527(1RAS 04360+2557) 


04:39:53.6 +26:03:05.5 


O-I 


WCCC, I ~ 85 


IRC-L1041-2 


20:37:21, +57:44:13 


O-I'' 


nearly edge-on'^ 


IRAS04302+2247 


04:33:16.45, +22:53:20.7 


O-I 


~ 90° 


ASR41 


03:28:51.291, +31:17:39.79 


II 


i ~ 80° = 


2MASSJ1628137-243139 


16:28:13.7, -24:31:39.00 


II 


~ 86° ' 


HVTau 


04:38:35.280, +26:10:39.88 


II 


multiple system (/ ~ 84° ^) 


HKTau 


04:31:50.900, +24:24:17.00 


II 


binary (i ~ 85° >■) 


UY Aur 


04:51:47.31, +30:47:13.9 


II 


~ 42° 



" Coordinate of the center of the 1' aperture. 
Inclination of the disk estimated by Tobin et al. j2008b . 

Judged from the similarity of its N3 image and spectrum to those of LI 527 and IRAS04302. 
Wolf, Padgett, & Stapelfeldt ( i2003i l 
= Hodappetal. (2004) 

' Grosso et al. ( 120031 1 

^ Inclination of HV Tau C (Monin & Bouvi erl2000t 
^ IncUnation of HK Tau B (Stapelfeldt et al. |T998] > 




1RAS04302 ASR41 





IRAS04302 



Fig. 1. (a) Images of our targets taken by AKARI with the N3 filter (2.7 - 3.8;um). The stellar light is dispersed from right (short 
wavelength) to left (long wavelength). The FOV is 1' x 1'. The coordinate at the center of the FOV is listed in Table[T] Images of 
IRC-L1041-2 and 2MASS J1628137-243139 are damaged by cosmic rays, while images of HV Tau, HK tau, ASR41, and UY Aur 
are affected by saturation and column pull down, (b) Spectral images of 1527 and IRAS04302. 



fourth image. Referring to this measurement, we selected images 
with a relative shift of within 0.5 pixels, in both spatial and spec- 
tral directions. Images with larger shifts were discarded. Cosmic 
rays were removed by subtracting the median of the possible 
combinations of three frames from the selected frames. We then 
stack and average these cosmic-ray-free spectral images. 

The local-sky region was selected by carefully avoiding neb- 
ulosity, field stars, and bad pixels. In the case of L1527, there is 
no 'sky' region, because the nebulosity extends all over the 1' 
X r aperture (Fig[T]). Since we performed slitless spectroscopy, 
the spectral image of the nebulosity is a convolution of the neb- 
ula spectrum with the brightness distribution, and should not be 
considered as 'sky'. We selected from the AKARI archive a 1' 



X r spectral image that can be used as sky; i.e. a dark region 
observed at a similar period and coordinate to those of L1527. 

After the sky subtraction, we defined the spectrum region of 
each star in the spectral image, referring to the coordinates of the 
objects and the N3 image. The signal was integrated over eight 
pixels (12 arcsec) in the spatial direction. As an alternative ap- 
proach. Noble et al. (in prep) fit the signal at each wavelength 
with the PSF and integrated the signal. We confirmed that the re- 
sult of the eight-pixel integration agrees well with this more so- 
phisticated approach. We note that the spectrum is not perfectly 
aligned with the rows of the pixels. We adopted the formulation 
as the tilt of the spectrum 

Y = -0.00746929(X - Xq) + Yq, (1) 
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where X and Y are the pixel number in the dispersion direc- 
tion and spatial direction respectively, and Yq is the center of the 
spectrum at an arbitrary position, Xq, in the dispersion direction. 

In Fig [U nebulosity is clearly seen towards LI 527, IRC- 
L 104 1-2, IRAS04302, and ASR 41. For these objects, we set 
io in Eq- (1) as the brightest position in the spatial direction at a 
wavelength (Xq) in the spectrum image, and integrated the spec- 
trum over eight pixels in the spatial direction at each wavelength 
(FigurelTJ. One exception is L1527; since the central point-like 
source is apparent and fainter than the nebulosity, we set as 
the central (faintest) position in the spatial direction, and inte- 
grated the signal over four pixels in the spatial direction (Figure 

m. 

The integrated signal was divided by the aperture correction 
factor and response function before being converted into a spec- 
trum. The response function was determined by the observation 
of a standard star (KF09T1). Wavelength calibration was done 
by the AKARI project team using the slit spectroscopy obser- 
vation of standard stars. In our observations, we determined the 
wavelength of the spectrum by measuring the shift between the 
slit and the direct light position of the target in the N3 image. 

According to the IRC team, the pixel-wavelength relation is a 
linear function. This linear correlation is confirmed in the wave- 
length range 2.5-4.5 fj.m, although the calibration is less certain 
at longer wavelengths. The CO absorption peak (~ 4.67 fim) is 
shifted to shorter wavelengths by 1-2 pixels in our data. Since 
similar shifts in this wavelength are reported in other IRC data, 
we shifted the data at > 4.5 yum by 1-2 pixels towards longer 
wavelengths. 

The error in flux was evaluated as follows. Firstly, we esti- 
mated the noise level in the 'sky' region and converted it into an 
error in flux using the response function. Secondly, we estimated 
the flux eiTors caused by the wavelength assignment; we calcu- 
lated fluxes by assuming that the shift between the slit and the 
direct light position of the target is 0.5 pixels larger/smaller. We 
then adopted the maximum of these errors. 

3.2. Absorption features 

To evaluate the ice absorption profile quantitatively, we normal- 
ized the spectrum with the continuum, which is assumed to be a 
second order polynomial, and was determined by referring to the 
three wavelength regions ~ 2.6 fim, 3.4 - 3.8 fim, and 4.2 - 4.9 
fim. The choice of wavelengths differed slightly among objects 
to ensure a good fit to all spectra. 

Template absorption profiles were obtained from "data 
base 2007" at Raymond & Beverly Sackler Laboratory for 
Astrophysicfl The database summarizes the absorption profiles 
of ice species measured in laboratories. Gerakines et al. (119951 
1 19961 1 measured the profiles and determined band strengths 
of various astrophysical ice analogs. From the database of 
Gerakines et al. (1995; 1996), we adopted profiles of pure ices 
at 10 K unless otherwise stated. In reality, interstellar ice is a 
mixture of various components such as H2O, CO, and CO2, and 
the peak position and shape of the absorption features change 
with the abundance ratio of the components. However, our spec- 
tral resolution is insufficiently high to perform such a detailed 
analysis of the ice mixture. 

The profiles of the 4.27 /urn CO2 band and the 4.67 fim CO 
band depend strongly on grain shape (Ehrenfreund et al l 19971 ). 
We adopted the continuously distributed ellipsoids (CDE) grain 
model for the CO band feature (Ehrenfreund et al .l997. Boogert, 

' |http://www.strw.leidenuniv.nl/~lab/| 



Blake & Tielens l2002l Pontoppidan et al. l2003l l. For die CO2 
band, on the other hand, we adopted the spectrum of Gerakines 
et al. (II995I 1, which is uncorrected for the grain shape, since it 
gives a good fit to our observational data. The pure CO2 fea- 
ture with grain-shape (CDE) coiTections agrees with neither our 
data nor the CO2 ice feature observed by ISO towards high-mass 
YSOs, which indicates that CO2 ice is not pure in interstellar 
medium. The ISO observations of the 4.27 jjm feature indicated 
that the CO2 ice is primarily in a polar (H20-rich) component 
of the ices (Whittet et al. 119981 ). However, to interpret our ob- 
servation, we also need to introduce an apolar component to re- 
produce the blue wing at ~ 4.2 /im. Since our spectral resolv- 
ing power is low, we simply use the spectrum of Gerakines et 
al. ( 119951 ) as a template, rather than combining multiple com- 
ponents of ice mixture with grain shape corrections. We also 
adopted the 3.05 fim H2O band of Gerakines et al ( 1995), since 
it is little affected by the grain shape. 

The spectral resolution of IRC, ~ 100, is not high enough 
to fully resolve narrow bands such as CO and CO2. We fitted 
these narrow bands using the following method: we considered 
the absorbance of ice from the database and convert it into a nor- 
malized flux spectrum, assuming the optical depth of ice along 
the line of sight. The spectrum was then convolved with the in- 
strumental line profile and re-binned to the resolving power of 
AKARI, to be compared with the observed spectrum. We vary 
the optical depth of ice until the best-fit is found by least squares 
fitting. Broad absorption features such as the H2O band at 3.05 
yum were not significantly altered by the convolution. 

In N3 images (2.7 - 3.8/zm), class O-I objects (L1527, IRC- 
L 104 1-2, and IRAS04302) are more extended than a point 
source. Among the three objects, LI 527 is the most extended, 
and IRAS04302 is the least extended. Since the IRC spectro- 
graph is slitless, the spectral resolving power depends sensitively 
on the spatial extent of the source. Assuming that the projected 
column density of ice in front of the extended sources is con- 
stant, the observed band profile is a convolution of the ice ab- 
sorptions with the spatial distribution of the light source in the 
dispersion direction. For LI 527, we integrated the central four 
pixels of the N3 image in the spatial direction around the source 
position to derive the spatial distribution of the light source in 
the dispersion direction (N3 profile, hereafter), with which the 
laboratory ice spectrum was convolved. This convolution better 
reproduces the observed features of CO and CO2 than the con- 
volution with the instrumental line profile. For IRC-L1041-2 and 
IRAS04302, however, we convolved the ice features with the in- 
strumental line profile, which reproduces the observed spectrum 
better than the N3 profile. Since the instrumental line profile at 
4 fim (N4 band) is naiTower than that of the N3 band, the N3 
profile might be too broad to be applied to CO and CO2 features. 
Unfortunately, imaging observation was not performed in the N4 
band. 

Once we determined the best-fit absorption feature, we used 
the original, pre-convolution profile to calculate the the column 
density of ice 

where A is the band strength (e.g. Whittet et al. 119931 ). We 
adopted the band strength listed in Table 2 of Gibb et al. (120041 ). 
unless otherwise stated (Table|2]). The CO ice column density of 
the CDE grain model was calculated to be 

A? = 6.03wA-', (3) 
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where A = 1 . 1 x 10"'^ cm molec ' and Tmax is the optical depth at 
the wavelength of the absorption peak (Pontoppidan et al. l2003l l 

Since the spectral resolution is low, and since some lines are 
saturated, it was difficult to estimate the eiTors in the ice column 
densities accurately. In the present paper, the column density er- 
ror was evaluated by simply calculating the error in the absorp- 
tion area. For CO2, we calculated the absorption area j rdv at 
4.2-4.3 i-im and its error caused by the flux error. Assuming that 
the absorption area had an error of 1 1 %, and that the CO2 col- 
umn density was estimated to be 10.0 x 10'^ cm"^, for example, 
the column density error was calculated as 1.1 x 10'^ cm"^. For 
H2O, we estimated the error by calculating the absorption area 
between 2.8-3.2 fim. If the red wing was much broader than the 
laboratory spectrum, we estimated the error from the absorption 
area at 2.8 - 3.0 yum. For L1527, towards which the H2O band is 
saturated, we used the wavelength range ~ 2.8-2.9 fim. The ab- 
sorption at ~ 4.7 yum is a combination of the absorptions of CO, 
XCN, and CO gas (§4.5). Here, the error is evaluated by calcu- 
lating the absorption area at 4.65-4.7 fim for CO and 4.57 -4.65 
fim for XCN. Uncertainties in ice column densities originating 
from the choice of the convolution function (instrumental line 
profile or N3 profile) are discussed in §4. 

The derivation of ice column densities described above im- 
plicitly assumes a simple geometry: a point light source behind 
a uniform absorbing material. Edge-on YSOs have more com- 
plicated structures. Stellar light scattered by the outflow cavity 
can be an extended light source (e.g. Tobin et al. 120081 ). even if 
it looks like a point source with a PSF of a few arc seconds, and 
the column density of the absorbing material (envelope and/or 
disk) also varies spatially. More accurate estimates of the ice col- 
umn densities requires the radiative transfer calculation of two- 
dimensional (2D) or three-dimensional (3D) models of envelope 
and disk structure (Pontoppidan et al. 2005). We compare the ob- 
served spectra with existing model calculations and discuss the 
location of ice in §5. 

4. Results 

4.1. Spectrum 

Figure |2] shows the derived spectra of our target YSOs. The spec- 
trum of L1041-2 resembles that of L1527, which suggests a sim- 
ilarity in the structure and evolutionary stage of these objects. 
The spectrum of IRAS04302, on the other hand, is less reddened 
than L1527. It might be in a more evolved stage than L1527 and 
L1041-2. 

2MASS J 1628 137-243 139 is very faint and is almost at the 
faintest limit for which we can derive a spectrum. The spectrum 
of ASR 41 shows a sharp decline at ~ 4.5 jum and is noisy at 
longer wavelengths because of the column pulldown caused by 
a bright star in the 10' x 10' aperture. 

HV Tau is a triplet system: HV Tau A and HV Tau B are a 
binary system with a separation of 0.0728" (Simon et al. 119981 ). 
and our target edge-on disk, HV Tau C, is separated by ~ 4" 
from the binary (Woitas & Leinert [19981 1. Since the FWHM of 
the AKARI PSF is ~ 4.3" in the NIR and the PSF is not circu- 
lar, it is difficult to disentangle the HV Tau C spectrum from that 
of the system. Hence, we integrated the spectral image over an 
eight-pixel width. The spectrum is dominated by HV Tau A and 
B at 2.5 - 4.5 fj.m, as they are brighter than HV Tau C by 4.22 
mag in the K band and > 1.71 magnitude in the L band. At 4.5 
fj.m, on the other hand, the contribution of HV Tau C may not 
be negligible, because HV Tau C is as bright as the HV Tau AB 
system in the N band (Woitas & Leinert il998i l. HK Tau is also a 




wave length [^m] wave length [lain] 



Fig. 2. The derived spectra of target YSOs in the NIR (closed cir- 
cles with error bars). The continuum (dashed line) is determined 
by the second-order polynomial fitting of the spectrum. 



binary system. The edge-on object, HK Tau B, is a faint binary. 
Unfortunately, it is impossible to disentangle the spectrum of HK 
Tau B from the primary, because the two stars almost align in the 
dispersion direction. It is disappointing that we were unable to 
discern the spectra of HV Tau C and HK Tau B. The absorption 
features are, however, still apparent in our spectra. We compare 
these absorptions with Subaru observations towards the HV Tau 
C and HK Tau B by Terada et al. ( 120071) in the following discus- 
sion. 

The dashed lines in Fig |2] depict the continuum. While the 
continuum varies between objects, reflecting their structure and 
evolutionary stage, we can clearly see absorption bands of H2O 
(3.05 /Km), CO2 (4.27 /im), and CO (4.67 fim) in most of the 
spectra. The assignment and fitting of each ice band are per- 
formed in the following subsections. 

4.2. H2O ice 

Figure |3] shows the spectra normalized by the continuum. The 
vertical axis is the optical depth of the absorption. Towards 
L1527 and IRC-L1041-2, the 3.05 fim water band is clearly sat- 
urated; the flux at the absorption peak is almost zero (Figure|2]i. 
We fitted the spectra in the range 2.9 - 3.1 fim with the labora- 
tory spectrum of pure H2O ice at lOK (Gerakines et al. |1995K red 
lines in Figure [3]). We also fitted the observed feature with the 
pure H2O ice feature at 160 K. The agreement is much worse 
than for the ice at 10 K, which suggests that the bulk of the ob- 
served H2O ice is not annealed. 
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Interestingly, the profiles fitted at 2.9 - 3.1 fim (solid lines in 
Figure |3]l reproduce the shallow hollow at 4.5 yum, which is the 
O-H combination mode (Boogert et al. 2000). While the 3.05 fim 
band is saturated, the 4.5 fim band suggests that the H2O column 
density cannot be significantly larger than we estimated. 

Table |2] summarizes the column density of H2O ice de- 
rived by fitting the 3.05 //m band. The class O-I objects L1527, 
IRAS04302, and IRC-L1041-2 clearly have much higher water 
column densities than the other targets, probably because of the 
larger mass of their circumstellar material and geometrical ef- 
fects (§5). 

The observed profiles show a red wing (> 3.1 yum), which 
is often observed in molecular clouds. The shape and depth of 
the red wing relative to the peak (3.05 fim) vary significantly 
among objects; it is especially deep towards L1527 and L1041-2. 
The wing is generally attributed to the scattering by large grains 
with ice mantles (e.g. Leger et al. |1983l l. The size of the grains 
responsible for the red wing is A/2n ~ 0.5 fim. 

ASR 41 has a double-peaked profile at the absorption peak 
~ 3 fim. We reduced each frame of the spectral image to confirm 
that the double-peaked profile is not caused by noise in a specific 
frame. The double-peaked profile remained, even if we changed 
the combination of frames in the stacking. We fitted the profile 
with (i) the H2O band alone (red line in Figure [3]) and (ii) the 
H2O band plus an NH3 band (2.98 fim) (blue line in Figure O. 
Although the blue line seems to fit the observed profile slightly 
better, it results in an unreasonably high column density of NH3 
(2 X 10^^ cm~^) compared with the H2O column density (7.8 x 
10^^ cm"^), and the bump at 3.1 fim is not fitted well. While 
the double-peaked profile seems to be robust, it could still be 
noise, since the flux at the absorption peak is only 1-2 mJy. We 
therefore prefer to fit the 3 fim absorption profile with the H2O 
band alone. 

In the polynomial fitting of the HV Tau continuum, we ig- 
nored the spectrum at > 4.5 fim, since the gradual rise of the 
spectrum at longer wavelengths may be caused by HV Tau C, 
while shorter wavelengths are dominated by HV Tau AB. The 
absorption band at 3 fim is much broader than the laboratory pro- 
file. While the red wing may be caused by scattering, the blue 
wing may be caused by the combination of the spectra of HV 
Tau A and B. The binary is close enough (~ 0.07") to be con- 
sidered as a point source in our observation, but the objects have 
similar brightnesses and different colors (Simon et al 1998). We 
tentatively fitted the absorption peak with the laboratory data to 
derive the water ice column density. 

Terada et al ( 120071 ) detected deep H2O absorption towards 
the edge-on disks of HV Tau C and HK Tau; they concluded that 
these absorptions originate in the disk ice. If there were no fore- 
ground clouds with ice, our spectra should be a combination of 
the edge-on disk with ice features and bright primaries without 
absorption. The observed 3 fim absorption of HV Tau is, how- 
ever, too deep to be accounted for by the disk ice alone; while 
the flux of HV tau C is 7.75 mJy at K band and 7.14 mJy at L' 
band, the 3 fim absorption is as deep as > 50 mJy. There must be 
a contribution of ice absorption in foreground clouds. The H2O 
column density in Table |2] should be considered as upper lim- 
its to the ice in foreground clouds, since we do not distinguish 
between the disk ice and foreground ice. However, the 3 fim ab- 
sorption of HK tau might be due to the disk ice of HK tau B; 
the flux of HK tau B is 15.9 mJy at 2.2 fim and 9.85 mJy at 3.8 
fim, while the absorption is < 20 mJy. If we were to assume that 
the flux of the primary is 200 mJy, subtract it from the spectrum, 
and fit the 3 fim absorption by the H2O ice feature, the ice col- 




wavelength [f.imj 



wavelength [|im] 



Fig. 3. Closed circles with error bars depict the YSO spectrum 
normalized by the continuum. The absorption feature at 3. 05 fim 
is fitted by the laboratory spectrum of pure H2O ice at 10 K (red 
line) (see text). 



umn density would be 3 X lO'^ cm ^, which is consistent with 
the value derived by Terada et al. (|2007] | (3.31 x 10'** cm"^). 



4.3. CH band at 3.5 fim 

While the long wavelength wing of H2O absorption is gener- 
ally attributed to scattering by large grains with ice mantles, it 
might also represent the stretching mode of C-H bond. Among 
various possible carriers of the C-H bond, CH3OH ice has long 
been observed via the substructures at 3.47 fim and 3.53 fim 
(e.g. Grim et al. [19911 Brooke, Seflgren & Smidi |T996] l. While 
the earlier works are mostly restricted to high-mass and inter- 
mediate YSOs, Pontoppidan et al (120031) identified the 3.5 fim 
feature towards three low-mass YSOs (out of 40 sources) as 
CH3OH ice. Boogert et al. (12008 1) detected the C-0 stretching 
mode of CH3OH at 9.7 fim towards 12 low-mass YSOs (out of 
41 sources) using Spitzer. Boogert et al. ( 1201 II) detected the 3.53 
fim band of CH3OH towards several field stars behind isolated 
dense cores; the CH3OH/H2O ratio ranges from 5 % to 20 %. 
The ratio of CH3OH to H2O column densities varies consider- 
ably among objects. 

Another possible carrier of the C-H bond is CH4, whose 
absorption bands at 3.32 fim and 7.67 fim have been detected 
towards high-mass YSOs (e.g. Lacy et al. 119911 Boogert et al. 
2004). Oberg et al. (120081) detected a 7.7 fim absorption feature, 
which is attributed to CH4, towards 25 sources out of 52 low- 
mass YSOs. Since CH3OH and CH4 are considered to be the 
precursors of large organic molecules in hot corinos and WCCC, 
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Table 2. Absorption peak, band strength, and column densities of ices towards objects 











cu 


ACJN 


Absorption peak [/im] 


3.05 


4.07 


4.27 


4.67 


4.62 


Band strength 10"''' cm molecule"' 


20 


4.3 


7.6 


see text 


5.0 


Column density 


10" [cnr^] 


10" [cm--] 


10" [cm-^] 


10" [cm--] 


10" [cm--] 


L1527(04368+2557) 


47±13 


0.96 (2.0%)" 


10±1.1 (21%) 


18+2.6 (38 %) 


1.5 +0.22 (3.2 %) 


IRC-L1041-2 


39±6.6 


3.7 (9.5%) 


9.4±1.6 (24 %) 


18+4.0 (46%) 


2.5+0.48 (6.4 %) 


IRAS04302+2247 


24±0.82 


5.3 (22%) 


6.8±0.64 (28 %) 


3.1+0.55 (13 %) 


0.53+0.064 (2.2 %) 


ASR41 


7.8±2.2 










2MASSJ1628137-243139 


6.7±2.6 










HVTau 


5.4+0.23 


1.0(19%) 


0.72 ±0.11 (13 %) 






HKTau 


2.1±0.23'' 




0.9±0.11 (43 %) 






UY Aur 


0.61±0.20 




0.54±0.086 (89 %) 







" Abundance relative to H2O ice. 
If the absorption originates in HK Tau B, ice column density is 3 x 10'* cm-^ (see text). 



respectively, it is important to constrain their abundances using 
the spectrum around 3.5 fim. 

Figure H shows the spectra of L1527, IRAS04302, and HK 
Tau around 3.5 yum. While the spectrum of L1527 is smooth, 
those of IRAS04302 and HK Tau contain absorption features. 
These features, however, do not match either CH3OH or CH4. 
A potential complication to the identification of, in particular, 
the CH4 3.32 jum band and to a lesser degree the CH3OH 3.53 
fim band, is the possible presence of the 3.3 fim PAH emission 
feature and the HI PfS line at 3.297 yum. We can indeed see an 
emission feature at 3.3 yum in HK Tau that could be ascribed to 
the PAH band. In the following, we derive upper limits to the 
column densities of CH3OH and CH4, keeping in mind that the 
low spectral resolving power of our spectra may blend adjacent 
emission and absorption bands. 

The thick solid lines in Figure |4] depict the upper limits on 
the CH4 absorption depth. We convolved the pure CH4 ice fea- 
ture at 10 K (Gerakines et al. 119951 1 with the instrumental line 
profile (IRAS04302 and HK Tau) or the N3 profile (for L1527). 
The upper limits to the CH4 column density are 1.6 x 10^' cm"^ 
(larger than the H2O column density), 3.1 x 10^^ cm'^ (13 % 
relative to H2O), and 1.6 x lO''^ cm-^ (76%) towards L1527, 
IRAS04302, and HK Tau, respectively. Owing to the low flux, 
the upper limit towards LI 527 is especially high. 

The dashed lines show the upper limits to the CH3OH ab- 
sorption depths, using a laboratory spectrum for pure CH3OH 
at lOK (Gerakines et al. 1995), convolved with the instrumental 
line profile. While the 3.53 fim feature is common among pure 
CH3OH ice and mixed ice (H20-HCH30H), the shape and depth 
of the 3.4 jum feature depend on the composition of the mixture 
(e.g. Pontoppidan et al. 2003 ). Hence, we derived upper limits 
to the CH3OH column density of the 3.53 //m feature. The upper 
Umits to the CH3OH column densities are 1.2 x 10'^ cm-^ (26% 
relative to H2O), 1.0 x 10"* cm-^ (42%), and 3.1 x 10'^ cm-^ 
(larger than the H2O column) towai'ds L1527, IRAS04302, and 
HK Tau, respectively. 



4.4. HDO ice 

In some objects, such as IRAS04302, we can see absorption 
around 4. 1 /im, where the OD stretching mode of HDO is ob- 
served in the laboratory (Dartois et al. 2003I I. Since the de- 
tection of HDO ice is very rare (Teixeira et al. 119991 1 and has 
not been confidently confirmed, we checked the response func- 
tion carefully and confirmed that the ~ 4.1 jjm feature is not 




30 

^25 
I20 



E 15 



10 

5 









IRAS04302 




3.2 3.3 3.4 3.5 
wave length [|im] 

Fig. 4. Three-micron spectra of LI 527, IRAS04302, and HK Tau 
(gray dots with error bars). The thin solid lines in LI 527 and 
IRAS04302 are the local (3.1-3.7 yum) continua fitted by second 
order polynomials, while the line in HK Tau is the best-fit ab- 
sorption of H2O (Figure |3]l. Solid and dashed lines depict the 
upper limits to the pure CH4 ice and pure CH3OH ice absorption 
depth, respectively. 
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Fig. 5. Spectra around the HDO absorption band. The soHd Hne 
is the Gaussian with FWHM of 0.2 jjm approximating the amor- 
phous HDO ice feature, while the dashed Hne approximates the 
annealed HDO feature. 



caused by an artifact in the response function. Figure|5]shows the 
spectra of this wavelength region towards L1527, IRC-L1041- 
2, IRAS04302, and HV Tau, which are fitted with a model of 
the amorphous HDO feature at 10 K: a Gaussian profile peak- 
ing at 4.07 yKm with a full-width half maximum (FWHM) of 0.2 
fj.m (solid lines). The spectrum of L1527 is not fitted well with 
this Gaussian; the absorption has a peak at longer wavelength 
(~ 4. 1 3 fim) and a narrower band width, which resembles an an- 
nealed, rather than an amorphous, HDO feature (Dartois et al. 



12003 ). We fitted the L1527 spectrum with a Gaussian peaking at 
4.13 fj.m and FWHM of 0.1 fim (dashed line in Figure|5]i. 

Although IRAS04302 has the deepest and smoothest absorp- 
tion around 4. 1 fim, this result should be taken with caution. The 
spectrum also shows a broad absorption at 4.5 - 4.6 fim, which 
cannot be fitted well by the absorptions of CO and XCN (see 
§4.6). These two absorptions (4.1 fim and 4.5 - 4.6 fim) could 
be related; there could be an alternative explanation, rather than 
HDO, CO, and XCN. The features themselves are, however, ro- 
bust. We have two independent data sets of IRAS04302, and the 
two broad absorptions appear in both data sets. 

We integrated the fitted spectra to derive the column den- 
sity of HDO (Table |2]l. The band strength was assumed to be 
4.3 X 10"'^ cm molecule"' (Dartois et al. 2003). Considering the 
small bumps and hollows that deviate from the HDO feature and 
the above discussion of IRAS04302 spectrum, the HDO column 
densities should be interpreted with caution. However, at face 
value, the HDO/H2O ratio ranges from 2 % (L1527) to 22 % 
(IRAS04302). Except in the case of L1527, the ratios are much 
higher than those obtained in the previous observations and the- 
oretical works: HDO /H2O < 3 % (Darto is eta l. l2003l Parise et 
al. 120031 Parise et al. l2^ Aikawa et al. lIOOST l. 

4.5. CO2 ice 

The CO2 band at 4.27 //m is clearly detected towards our tar- 
gets, except ASR41 and 2MASS J1628137. To subtract the H2O 
band at 4.5 fim, we divide the normalized spectrum by the best- 
fit H2O profile (red line in FigO. The absorption at 4.27 fim is 
then fitted by the laboratory CO2 profile convolved by the instru- 
mental line profile (N3 profile for L1527) (solid line in Fig|6l). 
To illustrate how the profile is broadened by the instrumental line 
profile, the laboratory profile before convolution is shown in the 



panel of IRAS04302 with the dashed line. While we adopted the 
pure CO2 profile measured by Gerakines et al. (1995) at 10 K, 
the actual shape of the CO2 profile depends on various factors, 
such as its mixing with other ice species, grain shape, and an- 
nealing (Ehrenfreund et al 1997, van Broekhuizen et al. 20061 1. 
We tried to fit the absorption bands with those of mixed ices 
of C02:C0=1:1 at 15 K, and H20:C02= 100:14 at 10 K, but 
the fits were worse than the pure CO2 ice of 10 K. Although 
the observed profile might be better fitted by other mixtures and 
temperatures, we do not pursue this fitting, because the spectral 
resolution of the IRC is not high enough. 

In regions close to the YSO, CO2 sublimates at ~ 70 K. We 
note that CO2 gas has absorption lines in almost the same wave- 
length range as CO2 ice. Although the observed feature might 
be a combination of CO2 ice and gaseous CO2, the low spectral 
resolving power makes it difficult to distinguish between gas and 
ice. We tried to fit the observed feature with CO2 gas of constant 
temperature (20K, 70K, and 150K), but the observed features are 
better fitted by the pure CO2 ice. 

We obtained the CO2 ice column densities listed in Table 
12] The errors listed in Table |2] were derived from the errors in 
the absorption area. For class O-I sources, the convolution func- 
tion (the instrumental line profile or the N3 profile) might add 
errors to the column densities. For IRC-L1041, if we were to 
adopt the N3 profile, the CO2 column density should be larger 
by ~ 30% than listed in Table |2j the wing region is reproduced 
better, but the agreement in the peak region is worse than in Fig 
|6] For IRAS04302, convolution with its N3 profile makes the 
absorption peak more flat and thus in closer agreement with the 
observation, while the CO2 column density is Uttle changed. 

We also detected CO2 absorption towards HV Tau and HK 
Tau. The absorption could, at least partially, originate in the cir- 
cumstellar disks. The H2O absorption towards HK Tau could 
originate solely in the disk of HK Tau B (§4.2). If so, the de- 
tected CO2 ice towards HK Tau should also be in the edge-on 
disk, since it is unlikely that the foreground clouds contain CO2 
but no H2O. The CO2 absorptions of HK Tau and HV Tau are, 
however, deeper than 10 mJy, while the fluxes of HV Tau C and 
HK Tau B seem to reach a maximum in the K band, and are 
smaller than 10 mJy at 3.8 fim and probably at longer wave- 
length (§4.2). Hence, we derived the CO2 ice column density 
by normalizing the spectrum with the continuum plotted in Fig. 
|2] as if the CO2 ice originates in the foreground clouds of the 
primary. 

The relation between column densities of CO2 ice and H2O 
ice is plotted in Figure |7] Our class O-I stars show similar 
CO2/H2O ratios to those observed by Pontoppidan et al. (I2008I I 
towards embedded low-mass YSOs using SST. Towards HV Tau, 
HK Tau, and UY Aur, on the other hand, the column densi- 
ties of H2O ice and CO2 ice are much lower than observed by 
Pontoppidan et al. (2008), assuming that the detected absorption 
features all originate in the foreground clouds. The CO2/H2O ra- 
tios toward these class II objects, however, agree with the range 
of ratios for low-mass YSOs from Pontoppidan et al. ( I2008I I. 

We note that the '^C02 band at 4.38 fim is apparent in the 
spectra of L1527 and IRC-L1041-2. The laboratory CO2 profile 
used to fit the 4.27 fim feature contains '-^002 with the terrestrial 
ratio I2c/i3c=89. In L1527, we fitted the 4.38 fim feature (thin 
solid line in Figure|6]l with the laboratory spectrum. The absorp- 
tion area (J rdv) of the thin solid line is about two times larger 

than that of the thick solid line, indicating that the '-C02/'^C02 
ratio could be lower than the terrestrial value. However, when 
considering the S/N, we found that the observed 4.38 fim feature 
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Fig. 6. Spectra fitted by the convolved laboratory spectrum of 
CO2 ice (solid line). The dashed lines depict the CO2 feature 
before convolution. In the panel of L1527, the 4.38 ;um feature 
is fitted by the convolved laboratory spectrum of '^COa (thin 
solid line). 



is consistent with the thick solid line in Figure |6l which fits the 
4.27 yum feature; the terrestrial ratio '^C/'"'C=89 are not signifi- 
cantly inconsistent with our observational data. 

4.6. COandXCN 

We found that CO ice has an absorption band at 4.67 jjm. Figure 
[8]shows the spectra around this wavelength towards L1527 and 
IRC-L1041-2, IRAS04302, HK Tau, and UY Aur 

Both L1527 and IRC-L1041-2 clearly have double-peaked 
profiles. The blue peak coincides with the XCN band, and the 
red peak matches the CO ice band. Although the double peak 
can be reproduced by combining XCN and CO ice, the observed 
features have an additional component at ~ 4.7 - 4.75 |/m. We 
assume that this corresponds to CO gas because, owing to the 
low sublimation temperature (~ 20 K), CO can be more easily 
sublimated and excited near the protostar than CO2. With the 
low spectral resolution of IRC, gaseous CO vibrational bands 
cannot be observed as lines but become a broad double-peaked 
(P branch and R branch) feature. It is thus difficult to derive the 
CO gas column density from our data. However, it is important 
to subtract the CO gas absorption, since the R branch of the CO 
gas feature coincides with the XCN feature. The distance be- 
tween the P branch and R branch increases with temperature. 
We fitted the absorption at 4.7 - 4.75 fim with a CO gas profile 
at 70 K that had been smoothed to the low spectral resolution 
(dashed lines in FigurelS). 

After the subtraction of CO gas absorption, we fitted the 
XCN feature, assuming a Gaussian peaking at 4.62 fim with 
FWHM of 29.1 cm^' (Gibb et al. (20041) (dotted lines in Figure 
[8]l, and the CO ice feature using the CDE grain model of 
Ehrenfreund et al. ( 1997) (thin solid lines in Figure|8]l. The dot- 
dashed line in Figure [8] depicts the laboratory spectrum before 
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Fig. 7. The relation between column densities of H2O ice and 
CO2 ice towards our targets (crosses) plotted over Figure 6 of 
Pontoppidan et al. (2008); the diamonds are low-mass stars ob- 
served by Pontoppidan et al. ( I2008I I. and the open triangles are 
the background stars from Whittet et al. (I2007I I and Knez et al. 
(120051) . 



convolution, and the thick solid lines in Figure [8] shows the sum 
of the three fitted components. 

The absorption feature around 4.67 fim is also observed to- 
wards HK Tau and UY Aur, but the feature does not peak at 4.67 
fim. We tentatively assigned the absorption to gaseous CO at 70 
K and 150 K for HK Tau and UY Aur, respectively (thin solid 
lines in Figure [8]). 

Table |2] lists the estimated column densities of CO ice. The 
correlation between the column densities of CO ice and H2O ice 
is plotted in Figure |9] over the correlation towards field stars of 
Whittet et al. (2007). The gray square depicts the ice column 
density ratio towards class I protostar Elias 29 (Boogert 200(^. 
Our targets show higher CO/H2O ratios than Elias 29, proba- 
bly because our targets are edge-on objects. IRAS04302 shows a 
lower ratio of CO/H2O than field stars, indicating that its circum- 
stellar material is warmer and/or less dense than typical molecu- 
lar cloud gas. 

Assuming that the XCN is OCN , we adopt the band 
strength of A(XCN)= 5 x 10"'^ cm mor', and estimate the 
OCN" column density to be (0.5 -2.5) x 10'^ cm"- towards class 
O-I objects (Table |2]l. The column density ratio of OCN /H2O 
ranges from 2.2 % (IRAS04302) to 6.4 % (L1041-2). Weintraub 
et al. (119941 ) detected the XCN feature towards a T Tauri star 
RNO 91 and derived an upper limit to the XCN/H2O ratio 
of < 9 %. More recent works, on the other hand, have de- 
rived much lower XCN/H2O ratio towards low-mass YSOs. van 
Broekhuizen et al. (120051 ) observed low-mass YSOs with VLT; 
the column density ratio of OCN /H2O is < 0.85 %. Although 
CO gas can contribute to the absorption at ~ 4.62 fim in our spec- 
tra, its contribution is subtracted by referring to the absorption at 
4.7 - 4.75 fim. In laboratory experiments, thermal or UV pro- 
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Fig. 8. The spectra around the CO ice absorption band fitted by 
a combination (thick solid line) of CO gas (dashed line), XCN 
feature (dotted line), and CO ice on the CDE grain model (thin 
solid line). The dot-dashed line in the panel of L1527 depicts the 
CO ice feature before the convolution with the instrumental line 
profile. 



cessing of N-bearing ice produces OCN" (Schutte & Greenberg 
I1997I I. Since our targets are edge-on objects, our lines of sight 
may preferentially trace the surface region of the disks (or the 
dense torus region of the envelope) where OCN" is produced 
by UV irradiation. The OCN" column densities towards IRC- 
L1041 and 1RAS04302, however, should be taken with caution, 
since there is residual absorption in the 4.5 - 4.6 yum region (Fig 
[8]). If there are other unknown carriers at these wavelengths, the 
estimated OCN" column densities are upper limits. 

Both LI 527 and 1RC-L1041-2 show weak absorption at ~ 
4.78 yum. Although it is close to the peak of '^CO ice feature, 
the observed peak is at a slightly shorter wavelength than that of 
i^CO. 



4.7. OCS 

IRC-L1041-2 has an absorption feature at ~ 4.9 /im. This coin- 
cides with the OCS ice band, which can be approximated by 
a Gaussian peaking at 4.91yL(m with a FWHM of 19.6 cm ' 
(Gib b et al 2004; Palumbo, Geballe & Tielens [T9971 Hudgins 
et al. I1993I I. Although the feature is at the edge of the observ- 
able wavelength range of the IRC and should be treated with 
caution, it does not disappear even if we change the combi- 
nation of frames to be stacked. We fitted the spectrum with a 
Gaussian (Figure [TOli. Assuming a band strength of 1.7 x 10"'** 
cm molecule ' (Hudgins et al. 1993), the OCS column density 
towards IRC-L1041-2 is 6.2 x lO'^ cm ^ If the abundance of 
H2O ice relative to hydrogen nuclei is 1 x 10""*, the hydrogen 
column density in the line of sight is 3.9 x 10^^ cm"^, and the 
relative abundance of OCS ice to hydrogen is 1.6 X 10 ^ which 
is much higher than the abundance of gaseous OCS, 1 x 10"^, 
observed towards TMC-1 (Irvine et al. J987i) . 
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Fig. 9. The relation between the column densities of H2O ice and 
CO ice towards our targets (square). The crosses are ice columns 
towards field stars observed by Whittet et al. (2007). The gray 
square depicts the ice column density towards the low-mass class 
1 protostar Ehas29 (Boogert et al. l20001 l. 
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Fig. 10. Spectrum of IRC-L1041-2 at 4.9 jum. The sohd line is 
a Gaussian with FWHM 19.6 cm ' approximating the OCS ab- 
sorption feature. 



5. Discussion 

5. 7 . Location of the Ice 

Young stellar objects have a complicated structure: a protostar, 
circumstellar disk, envelope, and outflow cavity. Furthermore, 
some objects might be accompanied by either ambient or fore- 
ground clouds. In our observation of edge-on YSOs, the light 
source is not necessarily the central star The NIR light of 3 - 5 
//m can originate from, for example, a wall of the outflow cavity 
that scatters the stellar light. The light source can thus be spa- 
tially extended relative to the projected spatial scale of ice abun- 
dance structures, even if it looks like a point source (e.g. Tobin 
et al. 2008; 2010), and the absorbing material in the envelope 
and/or disk also has non-uniform spatial distributions. Although 
the radiative transfer modeling of such an envelope and disk sys- 
tem is beyond the scope of the present work, it is important to 
discuss the location of ice in our targets. 



5.1.1. Class 0-1 objects 

Among the three class 0-1 objects, LI 527 is the most well- 
studied. Its bolometric temperature is less than 70 K, suggesting 
that the object is deeply embedded in dense gas (Furlan et al. 
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I2OO8I 1. Tobin et al. (I2008I I constructed an axisymmetric model 
of the envelope and disk with an outflow cavity referring to the 
SST images and SED of the object. According to the model, the 
"central object" we see in the N3 image is scattered light from 
the "inner cavity" with a size of ~ 100 AU. If so, the ice absorp- 
tion bands in our spectrum most likely originate in the envelope 
rather than the disk. The envelope density ranges from about 10^ 
cm^^ to 10^ cm^^ at ~ 100 AU to 15000 AU. The integrated en- 
velope column density in the line of sight is ~ 1 g cm"^, which 
corresponds to a hydrogen column density of A^h ~ 4 x 10^^ 
cm^^. This would, however, be an overestimate, because the en- 
velope structure is not axisymmetric in reality. 

Owing to the similarity of the continuum spectra, we expect 
L1041 and L1527 to have a similar physical structure. Varying 
depths of ice absorptions among these objects would thus probe 
the difference in the ice composition in the envelope. 

The SED modeling of the YSOs suggests that the envelope 
of IR AS04302 is sHghtly less dense than LI 527 (Furlan et al. 
I2OO8I 1. We found smaller column densities of ices and lower 
CO/CO2 ice ratios towai'ds IRAS04302 than L1527, which also 
indicates that IRAS04302 is more evolved and has a less massive 
envelope. 

Although our observations probe ice in envelopes rather than 
disks, the envelope material along our lines of sight (perpendic- 
ular to the rotation axis and outflow) is more likely to end up 
in the forming disks than be dissipated in the clouds. The high 
CO/CO2 ice ratio indicates that the majority of gases in the lines 
of sight are still dense and cold. The CO2/H2O ratio towards 
our objects is similar to the ratio observed by Pontoppidan et al. 
( 1200 8V who observed low-mass YSOs with various inclinations; 
the CO2/H2O ratio does not seem to vary with the inclination an- 
gle. 

The CO gas absorption and deep absorption feature of 
OCN , on the other hand, would originate in the warm (T > 20 
K) and irradiated region near the protostar However, other signs 
of heating by the protostar, such as crystalline feature of H2O ice 
(Schegerer & Wolf 120 10!) and CO2 gas, are not found in our data. 
We may need higher spatial resolution and spectral resolution to 
detect these features. 

5.1.2. Edge-on class II objects: ASR41 and 2MASS 
J1628137 

We estimated the H2O ice column density from our observations 
to be 7.8x10'^ cm"^ towards ASR41 and 6.7x10'^ cm"^ towards 
2MASS J1628137. These values are apparently much lower than 
we expect from the theoretical prediction of high H2O ice abun- 
dance n(H20)/nH ~ 10 in the disk midplane. The shallow ab- 
sorption is, however, caused by a geometric effect. Pontoppidan 
et al. ( I2OO5I 1 predicts that the 3 fj.m water band is the deepest for 
the disk inclination of ~ 70° and becomes shallower for higher 
inclination angles. If the inclination angle is > 72°, the optical 
depth of the band is predicted to be ~ 0.6, which is compara- 
ble to our observation of ASR41. Since the light source (scat- 
tered stellar light) is more extended than the absorbing material, 
Eq, (2) underestimates the ice column density towards edge-on 
disks. 

Possible contributions to the observed ice column densities 
from foreground clouds should also be considered. In the Taurus 
molecular cloud, the average H2O abundance is about 8.6 x 10"^ 
relative to hydrogen nuclei at Ay > 3 mag (Whittet et al. 119931 
Chair et al. |201 It . The H2O ice is not detected at lower Ay, which 
indicates that H2O ice is easily destroyed by photolysis or pho- 
todesorption near the cloud surface. Hodapp et al. ( 120041) esti- 



mated that the density and size of the molecular cloud around 
ASR 41 are 2.0 X 10"^° g cm"^ and 10"* AU, which corresponds 
to Ay ~ 1 mag. Although the relative abundance and the thresh- 
old Ay varies from cloud to cloud, the column density of the 
foreground gas towards ASR41, Ay ~ 1 mag, is too small to con- 
tribute significantly to the H2O ice absorption band. Therefore, 
the observed H2O ice most likely originates in the disk around 
ASR41. 

The foreground visual extinction towards 2MASSJ1628137, 
on the other hand, is estimated to be Ay = 2. 1 + 2.6 mag (Grosso 
et al. 2003 ). Assuming that the H2O ice abundance and thresh- 
old Ay are the same as in Taurus, the upper limit of the fore- 
ground H2O ice column density is ~ 2 x 10'^ cm"^, which is 
smaller than the observed H2O ice column (6.7 x 10'^ cm"^). 
Hence, we conclude that most or all of the H2O ice observed 
towards 2MASSJ1628137 originates in the disk. 

5.1.3. Other class II objects 

HV Tau and HK Tau are multiple systems. As the separation be- 
tween the multiple components in these systems is smaller than 
the size of our PSF in the NIR, we cannot extract the spectra of 
the faint edge-on objects. Terada et al. (2007) observed 2 - 2.5 
fim and 3-4 jjm spectra of these edge-on objects using Subaru 
and detected deep water absorption (t ~ 1 - 1 .5). The continuum 
flux of HV Tau C is ~ 7.75 mJy at 2.2 fim and ~ 7.14 mJy at 
3.8 fim, while the flux of HK Tau B is ~ 15.9 mJy at 2.2 yum and 
~ 9.85 mJy at 3.8 fim. Apparently, our spectra are dominated by 
the bright primaries. Both H2O and CO2 ice absorptions towards 
HV Tau should originate in both the disk and foreground clouds, 
since the depths of the absorption bands are larger than the con- 
tinuum flux of the HV Tau C. Towards HK Tau, the depth of 
the H2O absorption is comparable to the continuum flux of HK 
Tau B; if the H2O absorption originates in the disk around HK 
Tau B, the H2O ice column density is ~ 3 x 10'^ cm"^, which 
is comparable to the value obtained by Terada et al. (12007 ). The 
CO2 absorption, on the other hand, is deeper than the continuum 
flux of HK Tau B, and thus should originate at least partially in 
the foreground clouds. We also detected shallow H2O ice bands 
towards UY Aur, which again originate in the foreground com- 
ponent, since the inclination of the disk is ~ 42 degrees. 

If the H2O absorptions towards HV Tau and HK Tau orig- 
inated in the foreground clouds, the estimated H2O ice column 
densities would be smaller than the H2O ice column densities 
in the edge-on disks obtained by Terada et al. (2007) by about 
an order of magnitude; the H2O ice column density would be 
3.31x10'** cm-- for HK Tau B, and (2.69-4.20) x lO'** cm^^ for 
H V Tau C. Our observations thus support the argument of Terada 
et al (2007 ) that the disk component overwhelms the foreground 
component towards HK Tau B and HV Tau C. 



6. Conclusions 

We have observed ice absorption bands at 2.5 - 5 fim towards 
eight low-mass YSOs: three class O-I protostellar cores with 
edge-on geometry, two edge-on class II objects, two multiple 
systems with edge-on class II, and one not-edge-on class II ob- 
ject. 

Towards the class O-I objects, L1527, IRC-L1041-2, and 
IRAS04302, we have detected abundant H2O, CO2, and CO ice 
in the envelope. The column density ratio of CO2 to H2O ice is 
21 - 28 %, which coincides with the ratio observed by SST to- 
wards YSOs with various inclinations. The weak absorption at 
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~ 4.1 yum can be fitted by HDO ice; the HDO/H2O ratio ranges 
from 2 % to 22 %. The absorption in the vicinity of the CO 
band (4.76 yum) is double-peaked and fitted by combining CO 
ice, OCN , and CO gas. The large column density of CO ice sug- 
gests that the envelope is still very dense and cold, while OCN" 
and CO gas features would originate in the region close to the 
protostar The column density of OCN" is as high as 2 - 6 % rel- 
ative to H2O, which is much higher than previous observations. 
Our lines of sight (high inclinations from the rotation axis) may 
preferentially trace the regions with high UV irradiation, such as 
the surface of a forming disk and/or torus envelope. 

The spectrum of 1RAS04302 includes the 3.5 jum absorp- 
tion band, but the feature does not match either CH3OH or CH4. 
An OCS absorption band is tentatively detected towards IRC- 
L1041-2. 

Towards the edge-on class II objects, ASR41 and 2MASS 
J 1628 137-243 139, we have detected the H2O band. The low 
optical depth of the water feature is due to geometrical effects 
(Pontoppidan et al. I2005I) . The detected water ice mainly origi- 
nates in the disk. 

Both HK Tau B and HV Tau C are edge-on class II objects 
in multiple systems. In our spectra, which are dominated by the 
primaries, we have detected the absorption of H2O ice and CO2 
ice. Ices in both the edge-on disks and foreground clouds would 
contribute to the absorption, although H2O absorption towards 
HK tau could originate solely in the disk. Even if the observed 
features are due to the foreground clouds, the H2O ice column 
densities in the clouds are much smaller than those observed to- 
wards HK Tau B and HV Tau C with Subaru (Terada et al. l2007b . 
which confirms that the ice columns of the latter must originate 
in disks. The foreground H2O ice and CO2 ice are also detected 
towards UY Aur, which is not an edge-on system. We have ten- 
tatively detected CO gas towards HK Tau and UY Aur 
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